Unconventional pairing in the iron arsenide superconductors 
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We use magnetic long range order as a tool to probe the Cooper pair wave function in the iron 
arsenide superconductors. We show theoretically that antiferromagnetism and superconductivity 
can coexist in these materials only if Cooper pairs form an unconventional, sign-changing state. 
The observation of coexistence in Ba(Fei_ :E Co a ;)2As2 then demonstrates unconventional pairing 
in this material. The detailed agreement between theory and neutron diffraction experiments, in 
particular for the unusual behavior of the magnetic order below T c , demonstrates the robustness 
of our conclusions. Our findings strongly suggest that superconductivity is unconventional in all 
members of the iron arsenide family. 

PACS numbers: 74.70.Xa ; 74.20.Rp ; 74.25.Dw 



The determination of the mechanism of superconduc- 
tivity in the recently discovered[l|, Q iron arsenide com- 
pounds remains the prime issue in this field, requiring 
knowledge of the symmetry and of the internal structure 
of the Cooper pair wave function. A promising candidate 
for the pairing symmetry is the s" 1 -state ^proposed by 
electronic theories for superconductivity [3- 6], where the 
Cooper pair wave function changes sign between differ- 
ent sheets of the Fermi surface. In distinction to the case 
of (i-wave pairing^, H| in the much studied cuprates, no 
additional symmetry is broken for -pairing [9fl, mak- 
ing proposals for the determination of the wave function 
through interference experiments U 12 1 more complex 
and less conclusive. Important clues about superconduc- 
tivity in strongly correlated electron systems can, how- 
ever, be deduced by investigating their phase diagrams 
and the competition between different phases [13L Il4j. 

The iron arsenide superconductors manifest a rich 
phase diagram where antiferromagnetic (AFM), tetrago- 
nal (Tet), orthorhombic (Ort) and superconducting (SC) 
order are found in close proximity [ 1 5H2 1| . For some com- 
pounds, the transition between AFM and SC is of first 
order 15Wl7l| with regions of inhomogeneous phase coexis- 
tence. However, in Ba(Fei_ a; Co :E )2As2, experiments[18- 
[2fi | have established homogeneous coexistence of SC and 
AFM for intermediate ^-values. As the system enters 
the SC state, the ordered magnetic moment smoothly 
decreases with decreasing temperature[25l [2fj|. a behav- 
ior in sharp contrast to what is known for many conven- 
tional superconductors |27j, where AFM, associated with 
localized spins, can easily coexist with SC. This provides 
strong evidence for the fact that superconductivity and 
magnetic long-range order compete for the same elec- 
trons. 

In this Rapid Communication we demonstrate that the 
Cooper pair wave function in the iron arsenides is re- 
vealed via the coexistence and competition between su- 
perconductivity and magnetic order. We find that AFM 
and conventional phonon-mediated SC can not coexist, 



while unconventional -pairing is located near the bor- 
derline between phase coexistence and mutual exclusion. 
Therefore, the two phases can coexist only if Cooper pair- 
ing is unconventional with a sign-changing pairing wave 
function, whereas the absence of coexistence in other 
pnictide superconductors can not be used as evidence for 
conventional pairing. Our neutron scattering measure- 
ments confirm all aspects of our theory, including the 
novel re-entrance of the non-magnetically ordered phase. 
These findings strongly suggest that superconductivity is 
unconventional in all members of the iron arsenide fam- 
ily, given their similar electronic structure and transition 
temperatures. 

Microscopic model: we use a few basic ingredients to 
describe the main features of the iron arsenides: the elec- 
tronic structure is characterized by two sets of Fermi sur- 
face sheets, a hole pocket around the center of the Bril- 
louin zone and an electron pocket shifted by the ordering 
vector Q with Hamiltonian: 
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We use a circular hole Fermi surface, with £ p i = e^o — 
p 2 1 '(2m) — /i, and an elliptical electron Fermi surface, 
with £p+Q, 2 = -£2,o +pI/ (2m x ) +p 2 y / (2m y ) - /i. For the 
electron-electron interaction, H- m t, we include a magnetic 
electronic interaction I, i.e. i"Sp,p',q^ Ps i 'ss'V'p+qs'2 ■ 
ippi S 2 (T ss''4>p'~qs'i and a pairing interaction Vw, i.e. 

E P ,p', q ,H' v "»'V'p +qt ^ip i ;V'-p'-qii'V'p'ti' ; Although our 
key results are valid for arbitrary pairing matrix Vw , 
hereafter we will focus on the case of a predominant inter- 
band pairing Vw = V (1 — 6w). Depending on the choice 
for the sign of V, we consider the s H state, as arising 
from an electronic pairing mechanism with V > 
or the s ++ state that would result from electron-phonon 
interaction (V < 0). In the latter case, the Cooper pair 
wave function has the same sign in all Fermi sheets. We 
analyze the resulting model within a weak coupling mean 
field theory [11, [H and obtain the free energy density of a 
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Figure 1: Extrapolated zero temperature ordered moment 
M (T — 0, x) as function of doping x (panel a) and as function 
of Tn (panel b) for Ba(Fei_ a ;Co x )2As2. Points correspond to 
experimental data, whereas the solid line is the result of the 
calculation described in the text. In the insets, the red circle 
(blue ellipse) denotes the hole (electron) Fermi pocket. 



superconductor with antiferromagnetic long range order: 
I E log (4 cosh (^V) .(2) 



p,a— ± 



2k B T 



The SC order parameters "Fi and ^2 of the two bands 
and the staggered moment M are obtained by minimizing 
F (M, $ a ), where N is the system size and E pa are the 
positive eigenvalues of a state with AFM and SC order: 
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and 



^1 T "2 T c p,l "T- C p+Q.2 
4 (4,1 + A l) ( £ P +Q,2 + A ") as wel1 as S P = 

e P ,i£p+Q,2 - A| FM /2). Here A; = 
V^i and A A fm = IM refer to the SC and AFM sin- 
gle particle gaps. 

Neutron diffraction experiments: the neutron diffrac- 
tion measurements were performed on the HB1A diffrac- 
tometer at the High Flux Isotope Reactor at Oak Ridge 
National Laboratory on a series of Ba(Fei_ ;E Co :r )2As2 
single crystals using the same spectrometer configuration 
and data analysis methods described in (2||. The mag- 
netic integrated intensities were determined from rocking 
scans through the magnetic peak at Q tct = (|, |,3) as 
a function of temperature and put on an absolute basis 
using the known mass of the samples and the magnetic 
diffraction from the parent compound, BaFe2As2, mea- 
sured under identical conditions. 

The magnetic moment at zero temperature in the ab- 
sence of SC, M (T = 0, x), shown in Fig. 1, was deter- 
mined by extrapolating the measured order parameter 
M (T, x) above T c using a power law fit to the data. The 
ratios of the integrated intensities of the (^,5,1) and 
(2 >3>^) ma g ne tic reflections were monitored to ensure 
that there was no change in the moment direction as a 
function of temperature and composition x. No addi- 
tional reflections, e.g. incommensurate magnetic satel- 
lites, were observed, in agreement with other work (2l| . 



In Figs. 2a and b, we present a systematic study of the 
temperature and composition dependence of the mag- 
netic order by neutron diffraction. The key result is 
shown in Fig. 2b where we see that the magnetic or- 
der parameter, for superconducting samples, peaks at T c 
and decreases for T < T c . Indeed, for x = 0.059, our data 
show that there is a reentrance of the paramagnetic phase 
(magnetic long-range order is completely suppressed), in 
agreement with the predictions of our theory (see be- 
low). The opening of the superconducting gap removes 
states at the Fermi surface that otherwise contribute to 
the ordered moment, leading to a reduction of the or- 
dered moment below T c . 

Magnetic order in the absence of superconductivity: 
in Fig. 1, we demonstrate that our model of itinerant 
magnetism provides a description of the magnetically or- 
dered state that is consistent with the neutron diffrac- 
tion data. The calculation of M (T = 0, x) is done by 
setting A; = 0, but with magnetic order caused by an 
electron-electron interaction / ~ 0.95 eV, chosen to yield 
Tn = 140 K at x — 0. The other parameters used were 
ci.o = 0.095 eV, £ 2 ,o = 0.125 eV, to = 1.32TO c i octron , 
m x = 2m and m y — 0.3m, which yield an evolution of the 
Fermi surface with doping consistent with what is seen 
by angle resolved photoemission spectroscopy (ARPES) 
[30j. The suppression of AFM upon doping arises from 
the detuning of the two Fermi surface sheets; the hole 
sheet around the Brillouin zone center shrinks and the 
electron sheet grows upon electron doping (see insets of 
Fig. 1). The carrier density, x, is fixed under the assump- 
tion that each Co adds one electron. This analysis fixes 
all our parameters, except for the pairing interaction V. 
The latter is chosen to yield T c ~ 25K for z-values where 
AFM vanishes, yielding \V\ ~ 0.46 eV at x ~ 0.062. 

Competing order: without AFM order, many prop- 
erties of s ++ and s H pairing states are quite similar. 
This changes dramatically once AFM and SC compete 
for the same electrons. For s" 1 -pairing the excitation 
energies E p ± are fully gapped, whereas for the s ++ -state 
nodes occur once Aafm > A1A2, at momentum val- 
ues p n given by E Pn - = [3l|. This is true even if 
the nonmagnetic SC state is fully gapped. For the case 
where Ai = A2; nodes are located at e p „,i = £ p „+q,2 
(i.e. where Bragg scattering due to AFM is large, see 
Ref . (3l| ) . Note, however, that in general nodes are not 
guaranteed to emerge: for example, in the nested case 
£ p i = — Ep+Q^ with I > V, the s ++ -state remains fully 
gapped. To judge whether these AFM-induced nodes for 
,s ++ -pairing are relevant, one needs to analyze the free 
energy of Eq|2j which determines whether or not the two 
phases are allowed to coexist. 

In Fig. 2 we compare our theoretical results for the 
phase diagram and the temperature dependence of M 2 
with experiments. In Figs. 2c and d we show the cal- 
culated phase diagram and the behavior of M 2 for the 
s~* -state, using the parameters discussed above. The 
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Figure 2: a, The phase diagram of Ba(Fei_ 2 ,Co a: )2As2 deter- 
mined from neutron diffraction (solid symbols) as well as bulk 
thermodynamic and transport measurements (open symbols, 
Ref.[3). b, The AFM order parameter squared measured 
via neutron diffraction as function of temperature. The right- 
hand panel shows data in an expanded scale. Mo = 0.87/ib 
is the ordered moment at T = x — 0. c, and d, The phase 
diagram and theoretical ordered moment, obtained for an un- 
conventional s" 1 pairing state (inset). Note the SC and AFM 
coexistence region, in green. Panels e and f are analogous to 
c and d, but for conventional s ++ pairing (see inset). Hetero- 
geneous coexistence of AFM and SC regions at the first order 
transition occurs in panel (e) in a very narrow regime (dark 
red). 



phase coexistence and detailed temperature dependence 
of the ordered moment agree well with experiment, in- 
cluding the narrow doping regime with reentrance of the 
the AFM transition line. This is clearly different for the 
,s ++ -state: in Fig.2e and f we show that the two phases 
can not homogeneously coexist and are separated by a 
first order phase transition. Thus, -pairing and mag- 
netic order compete but coexist microscopically, whereas 
both phases are mutually exclusive in case of s ++ -pairing. 

This conclusion is robust and independent of specific 
details of the model, as follows from a Landau expansion 
of EqFZ]with respect to the order parameters. This expan- 
sion is performed near the multicritical point Tn = 2c j 
where both phase lines meet and where the decision 



about coexistence versus mutual exclusion takes place. 
Figs. 3a, b and c show the Feynman diagrams that are 
responsible for the order parameters coupling coefficients. 
The key diagram that is responsible for the different be- 
havior of s ++ and pairing in Fig. 2 is shown in Fig. 
3c, corresponding to a term M 2 (\E , 1 \E , 2 + *2*i) in tlie 
energy, which is sensitive to the relative phase, 0, be- 
tween 'J'i and V E , 2- Two partners of a Cooper pair in one 
band are coherently scattered into the other band where 
they recombine. AFM is essential, as it supplies the mo- 
mentum transfer Q needed for the scattering process. 
While the process contributes to the total energy for ei- 
ther pairing state, the phase of the pair wave-function 
determines the sign of this contribution, causing the sen- 
sitivity of the phase diagram with respect to the internal 
structure of the Cooper pair wave function. 

To illustrate the physical origin (and generality) of our 
results, we first discuss a simple limit that allows for ana- 
lytic treatment. Since the low energy electronic structure 
of the iron pnictides is nearly particle-hole symmetric, 
we consider Tn = T c and assume particle-hole symmetry, 
i.e. e p = e P) i = -£ p+ q,2, implying * = *i = e l9 * 2 
and I = \V\. In this limit, the Landau expansion of EqJ5] 
(relative to the nonmagnetic normal state) yields: 



F = 



M- 



#1 



M' 



•5(0) |*| M 2 



(4) 

which is highly symmetric in the two order parameters. 



F depends on the phase through 5 (0) = § (f + cc 
and on the two coefficients a = 21 — ^ ^ tanh (^f) 



and u 



r 

4NT 



sech 2 (If.) [T sinh 



e p ] /ep > 0. 



Therefore, the following results are completely indepen- 
dent of further details of the band structure dispersion, 
allowing us to draw general conclusions about the phase 
diagram for different microscopic pairing states. 

If 5 > 0, the two ordered states are separated by a 
first order transition, while homogeneous phase coexis- 
tence and second order transitions only occur if 5 < 0. 
The s ++ state, with g ++ = 5 (0 = 0) — u, is deep in the 
first order transition regime. Interestingly, for particle- 
hole symmetry, the state, with g^ = 5 (0 — tt) = 0, 

is at the border between regimes of coexistence and ex- 
clusion. The special symmetry in EqU at g — is di- 
rectly related to the emergent SO (6) symmetry that was 
found in electronic theories for s' 1 pairing [5] , as dis- 
cussed in Ref.j32|, suggesting that this result holds be- 
yond weak-coupling. This demonstrates that, for the 
pnictides, -SC is compatible with AFM, while mag- 
netic order enforces ne w g ap- nodes or strongly reduces 
the gap for ,s ++ -pairing[31], impeding phase coexistence. 
It is also interesting to note the connection of expression 
((3]) with the SO (5) model proposed for the cuprates 14j . 

Now, moving away from the special case of particle- 
hole symmetry, we note that the inclusion of an infinites- 
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Figure 3: a, b, and c, Feynman diagrams responsible for 
the coupling coefficients between the SC and AFM order pa- 
rameters. Gi (k) denotes the non-interacting single particle 
Green's function. Note that diagram c is sensitive to the rel- 
ative phase between and 4'2. d, Summary of the results 
for the coupling coefficient g, considering different band dis- 
persions. 



imal chemical potential /iora small ellipticity brings 

to small but positive values. However, when both the el- 
lipticity and ji are finite, g-\ can be negative. Using the 

parameters that lead to good agreement with the exper- 
imental results in Figs. 1 and 2, we find virtually the 
same result as before for s ++ , S±±. rj \ while for s" 1 
the coefficient assumes a negative value, ps —0.26, 
allowing for the phase coexistence presented in Fig. 1. 

We note that, for a d-wave state at particle-hole sym- 
metry, we obtain ^ = y| — ^, i.e. it is less compatible 

with AFM than the s H state. A similar result is ob- 
tained for the case of a nodal s H state. Even though our 
neutron diffraction measurements did not detect any in- 
commensurability, we checked that our main results still 
hold even for a small incommensurability. Furthermore, 
an extension of our model including the lattice degrees 
of freedom satisfactory describes the behavior of the or- 
thorhombic state below T c , as we show in [33j. 

The fact that the s H state is on the verge of co- 
existence and mutual exclusion with magnetism implies 
that the observation of different phase diagrams, with 
and without coexistence (l5l - l20j . does not imply different 
pairing states. There are situations where the s H state 
coexists with AFM, as in Ba(Fe 1 _. c Co. c ) 2 As2, and oth- 
ers where it does not - as, for example, in the case of 
two detuned circular bands [29] (see Fig. 3d). Thus, the 
presence of a first order transition does not imply s ++ - 



pairing, whereas the inverse is true: observing phase co- 
existence in the iron arsenides disallows s ++ SC. This 
makes Ba(Fei_ x Coa;) 2 As2 a crucially important member 
of the pnictide family. 

We thank Alfred Kracher for performing the WDS 
measurements and Sergey Bud'ko for assistance with 
the thermodynamic and transport measurements. This 
work was supported by the U.S. DOE, Office of BES, 
DMSE. Ames Laboratory is operated for the U.S. DOE 
by Iowa State University under Contract No. DE-AC02- 
07CHI1358. 



[9 
[10 

[11 
[12 

[13 



[14 
[15 
[16 
[17 
[18 
[19 

[20 
[21 
[22 

[23 
[24 
[25 
[26 

[27 

[28 
[29 

[30 
[31 



Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008). 
M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008). 

I. I. Mazin, D. J. Singh, M. D. Johannes and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 
K. Kuroki et al, Phys. Rev. Lett. 101, 087004 (2008). 
A.V. Chubukov, D. V. Efremov, and 1. Eremin, Phys. 
Rev. B 78, 134512 (2008). 

I. I. Mazin and J. Schmalian, Physica C 469, 614 (2009). 
D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 

C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

V. Stanev, J. Kang, and Z. Tesanovic Phys. Rev. B 78, 
184509 (2008). 

D. Parker and I. 1. Mazin, Phys. Rev. Lett. 102, 227007 

(2009) . 

J. Wu and P. Phillips, Phys. Rev. B 79, 092502 (2009). 
J. Linder, I. B. Sperstad, and A. Sudbo, Phys. Rev. B 
80, 020503(R) (2009). 

Y. Zhang, E. Dernier, and S. Sachdev, Phys. Rev. B 66, 
094501 (2002); S. Sachdev, Phys. Status Solid B 247, 537 

(2010) . 

S.-C. Zhang, Science 275, 1089 (1997). 
H. Luetkens et al, Nature Mater. 8, 305 (2009). 
J. T. Park et al, Phys. Rev. Lett. 102, 117006 (2009). 
T. Goko et al, Phys. Rev. B 80, 024508 (2009). 
N. Ni et al, Phys. Rev. B 78, 214515 (2008). 
J.-H. Chu, J. G. Analytis, C. Kucharczyk, and I. R. 
Fisher, Phys. Rev. B 79, 014506 (2009). 
F. Ning et al, J. Phys. Soc. Jpn. 78 013711 (2009). 
C. Lester et al, Phys. Rev. B 79, 144523 (2009). 
Y. Laplace, J. Bobroff, F. Rullier-Albenque, D. Colson, 
and A. Forget, Phys. Rev. B 80, 140501 (R) (2009). 
M. -H. Julien et al, Europhys. Lett. 87 37001 (2009). 

C. Bernhard et al, New J. Phys. 11, 055050 (2009). 

D. K. Pratt et al, Phys. Rev. Lett. 103, 087001 (2009). 
A. D. Christianson et al, Phys. Rev. Lett. 103, 087002 
(2009). 

W. Baltensperger and S. Strassler, Phys. Kondens. 
Mater. 1, 20 (1963). 

M. Kato and K. Machida, Phys. Rev. B 37, 1510 (1988). 
A. B. Vorontsov, M. G. Vavilov, and A. V. Chubukov, 
Phys. Rev. B 79, 060508(R) (2009). 

C. Liu et al, larXiv:0910.1799l f2009k 

D. Parker, M. G. Vavilov, A. V. Chubukov, and 1. I. 
Mazin, Phys. Rev. B 80, 100508(R) (2009). 

[32] D. Podolsky, H.-Y. Kee, and Y. B. Kim, Europhys. Lett. 



■5 



88, 17004 (2009). 
[33] S. Nandi el al, Phys. Rev. Lett. 104, 057006 (2010). 



